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Abstract 

The results obtained by Pauli in his 1926 article on the hydrogen atom made essential use 
of the quantum dynamical so(4) symmetry of the bound states. Pauli used this symmetry 
to compute the perturbed energy levels of an hydrogen atom in a uniform electric field 
(Stark effect) and in uniform electric and magnetic fields. Although the Stark effect 
on hydrogen has been studied experimentally, Pauli's results in mixed fields have been 
studied only for Rydberg states of rubidium atoms in crossed fields and for lithium atoms 
in parallel fields. 
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1 Introduction 

In 1926 there appeared nearly simultaneously two articles by Schrodinger |^ and Pauli 
that solved one and the same central problem of the newly discovered quantum mechanics: 
the computation of the hydrogen atom spectrum. In most textbooks on quantum mechanics 
the eigenvalue problem is solved using Schrodinger's approach with spherical coordinates and 
simultaneously diagonalizing the set of operators Hq, L ^, L^. The main benefit of this approach 
is that one can deal, in principle, with any central potential and that the eigenstates have 
definite parity, which is of practical importance in the study of the selection rules for radiative 
transitions. However, a simple problem like the Stark effect, is quite cumbersome to deal with 
perturbatively, except for the low lying states. 

An essential aspect of the Coulomb problem is that the Schrodinger equation is also sepa- 
rable in parabolic coordinates This super-separability is related to the conservation of an 
additional vector m, the so-called Runge-Lenz (RL) vector, which is specific to the Coulomb 
potential. Despite its name, the discovery of this conserved vector in Newtonian mechanics 
goes back to Hermann in the 18th century, as noted by Goldstein |p, who mentions also early 
works by Laplace and Hamilton. 

Pauli succeeded in defining a quantum extension M of the RL vector which is an observable, 
that is, an hermitian operator, and commutes with the Hamiltonian Hq. This success enabled 
him to calculate the spectrum of the hydrogen atom by an abstract approach. The quantum 
extension of the RL vector was used to obtain the dynamical so(4) symmetry of the bound 
states. This symmetry, combined with the general theory of angular momentum allows to 
express the energies of the hydrogen atom 0. 

The main drawbacks of Pauli's approach are that it works only for a Coulomb potential and 
that the eigenstates have, in general, no definite parity. However, in the same article, Pauli 
gave two non-trivial applications of his ideas: the first order perturbation of the energy levels 
under a uniform electric field (Stark effect) and in the more complicated case of uniform electric 
and magnetic fields. Guided by an analogous relation from classical mechanics, Pauli guessed 
the following crucial relation between quantum expectation values : 

{2HoR) = l{M), (1) 

where R is the position operator and the expectation value should be taken in any eigenstate 
of the Hamiltonian Hq. It is natural to call Eq. (|1]) Pauli's identity and we will later discuss 
its origin. 

Textbooks that discuss more advanced topics related to the quantum RL vector (see Refs. 
do not treat these interesting applications. The main aim of this note is to popularize 
Pauli's results. In Sec. II we introduce the quantum RL vector and use it to derive the spectrum 



of the hydrogen atom [|I0| . In Sec. HI we present a new derivation of the crucial identity 



2HoR = ^M+^[Ho,f], (2) 

which implies Eq. (p. It is interesting to note that Eq. (0) was first derived by a clever 
calculation of commutators by Becker and Bleuler [jll| fifty years after Pauli's article. Our 



derivation will make apparent that it is the natural quantum extension of a classical relation. 
In Sec. IV we consider the hydrogen atom in the presence of uniform electric and magnetic 
fields. We then use Pauli's elegant approach to obtain the first-order perturbed levels. We 
conclude with a short account of the experimental checks, some of which are quite recent. 
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2 Background material 

To set our notation (we stick closely to those of 0), let us recall the basic properties of the 
quantum RL vector. As far as possible, we will use upper case letters for quantum operators. 
We will write the hamiltonian of the hydrogen atom 

where /i is the reduced mass, q is the proton charge and k = g^/47reo. Although going from the 
classical angular momentum I = -r A p to the operator L = R A P is not ambiguous and leads 
to an hermitian operator L, the situation is somewhat more complicated for the RL vector. 
Classically its definition 

_, 1 ^ -> r 

rh = —pAl — K- (4) 
/i r 

shows that there are quantization ambiguities due to the lack of commutativity of the operators 
P and L. There are two possible (but non-hermitian) corresponding operators, PAL and 
(P A Ly = —L A P. Pauli observed that the simplest choice for M so that it is hermitian is 

M = —{P AL- LAP) - K^. (5) 
2/i R 

From these definitions we can check that both L and M are conserved at the quantum level, 
that is, that they commute with Hq. 

The classical Poisson brackets involving / and fh generalize to commutators of L and M: 



[Li, Lj] — ihtijk Lki 

(6) 



[Li,Mj\ = iheijkMk 



[Mi,M,] = (^-^^ -iheijkLk, 

Two further important relations should be noted: 

L- M = M ■ L = 0, 

The checks of the conservation of the RL vector and of Eqs. (|^) and (|^) are quite involved: a 
detailed calculation may be found in Ref. [§| (p. 462) and in Ref. (p. 265). 

Let us restrict ourselves to the subspace of the Hilbert space spanned by the bound states 
lip) with energy E <0. H^^ and {-Hq)-^/'^ can be defined by 

We also define 



1/2 

2Hn 



M'={—^\ M. 
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Because L and M commute with Hq, the algebra in Eq. becomes 

[Li, Lj] = iheijk Lk, 



(9) 



which we recognize as an so (4) algebra with generators L, M' . For further use we define the 
generators 

I = \{L + M'), 

. . (10) 
K=\{L-M'). 

Their commutation rules are quite instructive : 



[Ki, Kj] = iheijkKk, 



because they show that the operators / generate an so(3)/ Lie algebra, while the operators K 
generate an so{3)k Lie algebra. Furthermore these two algebras are completely independent 
(commuting), and their Casimir operators (i. e. operators commuting with the whole algebra) 
are for so(3)/ and K'^ for so(3)x- 

The correspondence given in ([To|) between L, M' and J, K gives just the change of basis 
that displays the isomorphism so(4) ~ so{3)j © so{3)k- From the general theory of quantum 
angular momentum 0, we can find a basis common to the operators J^, K^, J^, such 
that 

/ ^ \i, rrii; k, rrik) = i{i + 1)^^ \i, mf, k, mk) 



Iz\i,mi;k,mk) = rui h\i,mi; k,mk) 
\i, mf, k, iTLk) = k{k + l)h'^ \i, mf, k, m^) 



(12) 



mi G [—i, . . . , and G [—k 



with i,k = 0, |, 1, 
that the operators / and K are constrained by 

L- M' = 

J 2 = 



K^\i,mf k,mk) = mkh\i,mf, k,mk), 

, +k\. The first relation in (|^) implies 



(13) 



2{P + K'^) + h^ 



(14) 



while the second relation can be written as 

showing that the states \i,mf,i,mk), are also hamiltonian eigenstates. At the level of the 
eigenvalues, the operator relation in Eq. (0) gives 



1 2 1 

2^" ^ 



4^(z + i)h' + h' = {2i + lyn'. 



(15) 
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so that the inverse energies are given by the eigenvalues of the Casimir operators of so(4). Up 
to the identification of the principal quantum number n = 2i + 1, we obtain for the energy 
eigenvalues 

where n = 2i + 1 = 1, 2, . . . The eigenstates \i, rrii; i, m^) are such that rrii G [—i, . . . + i] and 
nik G [—i, . . . + i]. The degeneracy is easily seen to be (2i + 1)^ = n^. 

These states are not eigenstates of L^. Indeed the relation L = I + K shows that / may 
take any value between and 2z = n — 1, in agreement with the Schrodinger approach. Note 
that the complete set of commuting observables diagonalized here are Hq, I^, and Kz while in 
Schrodinger's approach one takes Hq, L^, and L^. Note also that the angular momentum L is 
axial, while M is a true vector, so that J, K have no definite parity, and hence their eigenstates 
cannot have any definite parity. This absence of a definite parity is also obvious from the fact 
that the states \i,mi;i,mk) have no definite value for the orbital angular momentum /. Let us 
conclude that even for parity invariant potentials the eigenstates need not be of definite parity 
as can be seen in the discussion of the simpler case of a particle in a box, given in Ref. [0 
(p.328). 



3 From classical to quantum Pauli's identity 

Equation (|I]) has a classical content, because it involves mean values. Indeed, for any classical 
periodic trajectory, the mean value of a dynamical quantity /(r), is defined by 

where T is the period and A = fir'^d(j)/dt is the angular momentum. For the Coulomb case, we 
have 

— = 1 + e cos 0, 
r 



1 + 2^, (18) 



T = 2n- 



A(l-e2)3/2' 

where tq = X^/fin, £ is the total energy and e the eccentricity of the orbit. We can show 
that the classical RL vector, lying in the trajectory plane, is directed from the origin to the 
perihelion, which lies on the x-axis. Its components (see Ref. p (p. 460) for the details) are 
given by 

rrix = ne, rUy = rriz = 0. (19) 

TIf we take / = r , we have 



{x) 



/iTq f^"^ COS0 SvTyuerQ 



AT (l + ecos0)3 Ar(l-e2)5/2' 

(20) 



urn f^"" sine/) 
AT /_„ 1 + ecosdV 
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Some algebra leads to 

2S{x) = ^,^e, {y) = {z) = 0. (21) 

If we compare Eqs. ([l9|) and (21) and identify the total energy with the classical hamiltonian 
ho, we obtain Pauli's relation in classical mechanics 

2/io (r=) = ^m. (22) 

Equation ( ^2]) is an equality between mean values of classical physical quantities, because 
the RL vector is conserved. It follows that the difference between these two classical physical 
quantities must be a total derivative (whose mean- value necessarily vanishes). So we need to 
find some vector r such that 

2ft„,= 5«_|f. (23) 

Since on dimensional grounds the sought vector has dimension of momentum x (length)^, it 
must have the structure 

T = a {r ■ p) r + br'^ p. (24) 
The coefficients a and b, fixed by enforcing relation (p3|), give for r the expression 

T = - {r ■ p) r — f'^ p, (25) 
which can be written more suggestively, using a Poisson bracket as 

The above relations involve classical (commuting) quantities ; to generalize Eq. ( pB]) to the 
quantum level, we try the following quantum extension 

2HoR = ^M+^[Ho,f], (27) 

where the quantities involved are now operators and T is the unknown quantum extension of 
the classical quantity r. We point out that in the left hand side of Eq. (p7\), the position of the 
hamiltonian operator with respect to R is important. In this case T can be written in terms of 
three vector operators 

f = aR{R- P) + b{R- R)P + ihcR. (28) 
If we require Eq. ( pTj) to be satisfied, we obtain the result |13 



f =^R{R- P) - {R- R)P + ihR. (29) 

Pauli's relation (|I|) follows from Eq. ( pTj) when we take its expectation value between eigenstates 
of Hq. We can write 

{2HoR) = 2E^^^{R) = ^{M). (30) 
For further use we will write Eq. (^) as 

R-1{2E(:Y'M, (31) 

The symbol ~ indicates that the equality holds only when expectation values between eigen- 
states (with energy En^ ) of the unperturbed hamiltonian are calculated. 
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4 Electric and magnetic fields 



In the approximation of an infinitely heavy nucleus, /i = m the electron mass. If we also neglect 
the diamagnetic terms, the perturbation due to the uniform electric field E and magnetic field 

— * 

B writes 



Hi = qE-R + —B-L, 
2m 



(32) 



where q is the proton charge. For the first order computation we need the matrix elements 
(i, 'm^■, i, m'^|ifi|i, mf, i, m^) for states of definie energy En'^ where i = {n — l)/2 has some fixed 
value. Its first piece, using Pauli's identity, in the form given by Eq. (|3T|), we have 



E-Rr^^{2E^^Y'E-M 



-2mEj^yr'^^E-M' = ~n^E-M', 



(33) 



with qq = h? / fiK. It is then sufficient to use 



M' = 1 -K, 



L = I + K, 



(34) 



to obtain 



H^r^ {--nS + B 



I 

h 



(35) 



with £ = qaQE,B = ^bB, and fj,B = qh/2m. Notice that dimensional analysis shows that £ 
and B have the same dimensions as energy. 
Let us now define 



±ln£ + B 



±ln£ + B 



E, 



±-n£ + B 



n"^ £^±3n£-B + B^, (36) 



so that the perturbation can be written as 

E. 



E^ 



H^^ — (u_ .l) + ^^V+.K). 



(37) 



This perturbation is made up of two completely independent (commuting) pieces which can be 
diagonalized separately. Indeed, one can check the relation |]T^ 



V 



(3J 



with 



z7 = sin 6{cos(j)i + sin j) + cos 6 k, 



u Ak 
\uAk\\ 



sm(p% — COS0J, 



which states that a rotation of angle 6 around the axis /2 transforms the vector u into the unit 
vector k along the z axis. It can be checked using the well-known identity 



Q ^-iA = 5 + , [^^ 5] + _ 5]] _^ 



01 
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and the so{3)i algebra for the generators /. It follows that the operator u ■ I has the same 
eigenvalues as 1^, i.e. rriih. Similar results are also valid for K. If we adapt this result to the 
perturbation given by Eq. (|37|) we can write 



u+-K = ei'^+^+-^ K, e-*^+^+-^, cos^+ = u+-k, /2+ 



v_ ^k 



\u- A k\ 

1^+ Ak 
\u+ A k\ 



(39) 



and we obtain the first order correction to the energies 

A E^^^ =mi- E_+mk- E+. (40) 

Recall that the quantum numbers m, and take all integer or half-integer values between 
-(n-l)/2and+(n-l)/2. 



4.1 The first order Stark effect 

For vanishing magnetic field, we can take the electric field along the z axis. From Eq. (|35D we 
obtain the first order perturbed energies 

AEW = -^gao 1^1 n(m,-mfc). (41) 

The levels are split into 2n — l sub-levels, each with residual degeneracy n— |m|. Eq. (^) is the 
celebrated first-order Stark effect formula, which was obtained in 1926 simultaneously by Waller 
ISf , Wentzel |l^, and Epstein by perturbation calculations in parabolic coordinates. A 



more handy reference is Ref. [Q, where perturbative results up to second order. The reader is 
urged to compare these extensive calculations (even in first order) with the elegance of Pauli's. 
The experimental results from Ref. |jl8| are given in figure 1 for the hydrogen atom. The 



low lying states exhibit a nearly linear field dependence up to a value of 10~^ au (atomic unit), 
but this value decreases with increasing n. This hnear regime accounts also for some crossings 
of the energy levels corresponding to different values of n. The second order corrections are 
also sizeable in this region for some levels. For higher field values the situation is much more 
complicated since the states are broadened by tunnelling effects or even ionized, and cannot be 
accounted for by perturbation computations around the bound states! This complicated region 
lies to the right of the solid curve. 



4.2 First order Zeeman effect 

For vanishing electric field, we can take the magnetic field along the z axis and we are back to 
the Zeeman energies 

AE(i) =fiB\B\ m, (42) 

where m = -|- m^. As for the Stark effect the energy levels are split up into 2n — 1 levels, 
each with residual degeneracy n — \m\. 
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4.3 Parallel versus crossed fields 

Pauli's result simplifies for the special case of crossed electric and magnetic fields, i. e. such 
that E ■ B = 0. The perturbed energies become 



A E^^^ = E^-m, with E± = ^^n^ {qao Ef + (/is B) \ (43) 

In some sense the system is still in a Zeeman-like regime, with the same splitting pattern, but 
with a more complicated mixed field dependence appearing in the factor E^^. Equation (^3]) has 
been checked for Rydberg states of rubidium atoms (large principal quantum number n = 34), 



which are essentially hydrogen- like |T^. The experimental results, taken from Ref. [|T9|, are 
reproduced in Fig. 2. The different lines correspond to different values of m. Note how small 
the window is for the electric field (between and 20 V/cm) and for the magnetic field (between 
and 6 10~^ T) to observe Pauli's quantization, but the results are in good agreement with 
( |45| ) for weak fields. 

An experiment dealing with parallel electric and magnetic fields has been described in Ref. 
^ for lithium atoms. These results are quite interesting, since the symmetry breaking is 



completely different from the crossed case, with no left over degeneracy. A quite large set of 
references to theoretical as well as experimental work may be found in [^]. 

The second order perturbative computation was done quite recently ||22[,[^. The most 
interesting aspects cover now the non-perturbative transition to chaos {B > lOT) which is 
quite harder to deal with. The interested reader should consult Ref. 
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FIGURE CAPTIONS : 



Figure 1: Stark effect for the hydrogen atom. An atomic unit (au) corresponds to an electric 
field of 5, 14 10^ v/cm. From Ref. 



IS 



Figure 2: Rydberg states of rubidium atoms in crossed electric and magnetic fields. The electric 
field unit is 1 V/cm and the magnetic field unit 10^^ T. From Ref. |]19| 
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